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Studying the jellium model in the Hartree-Fock approximation, Overhauser has shown that spin 
density waves (SDW) can lower the energy of the Fermi gas, but it is still unknown if these SDW 
are actually relevant for the phase diagram. In this paper, we give a more complete description of 
SDW states. We show that a modification of the Overhauser ansatz explains the behavior of the 
jellium at high density compatible with previous Hartree-Fock simulations. 
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The simplest model of electronic structure is jellium 
- electrons embedded in a homogeneous background of 
opposite charge such that the system is neutral. This 
model is a good starting point to describe properties of 
simple metals such as sodium[T]. At zero temperature, 
the only parameter of this model is the density n, or 
the dimensionless parameter = 3/(47rna^)^/^, where 
qb = fi^/{me^) is the Bohr radius. Within the Hartree- 
Fock approximation (HF), Overhauser has shown that 
the Fermi gas is unstable under a spin density wave 
(SDW) m [5] . Only recently, almost 50 years after Over- 
hauser’s prediction, explicit numerical estimates of the 
HF ground state have shown SDW evidences of the elec¬ 
tron gas in three [IHZ] and two dimensions [H |8]. Still, a 
quantitative estimate of the variation of the SDW am¬ 
plitude and energy in the high density region < 1 
is missing[S]. In this letter, we generalize Overhauser’s 
ansatz and provide a quantitative solution of this long¬ 
standing problem. 

The key point is to search for a solution in a non- 
perturbative way. Indeed, small domains exist around 
the Fermi surface where the one-body states differ rad¬ 
ically from a single plane wave. These states of wave 
vector k are coupled with the wave vector k -|- Qk where 
Qk is constant over each domain. The larger is this do¬ 
main, the larger will be the energy gain of the SDWs. 
One way to enlarge this domain is to cut the top of 
sphere as explained by Overhauser [5]. In the following 
we show that adding a small cylinder on the top of the 
truncated sphere as shown in Figcan increase the en¬ 
ergy gain of the SDW by orders of magnitude compared 
to Overhauser’s ansatz. Furthermore, we provide an ex¬ 
plicit estimate of the energy gain of the SDW state. As 
we will see, the optimal size of these domains dramati¬ 
cally shrinks with increasing density, resulting in a ex¬ 
tremely rapid decrease of the tiny SDW energy gain and 
explaining the difficulties of observing SDW in the high 
density region. 

Our semi-analytical results presented here are com¬ 
pared to recent HF results [6] obtained with periodic mod¬ 
els. Indeed, the Overhauser’s ansatz is in fact a periodic 
model (a crystal where the one-body states are limited 


to the first mode) as soon as the set vectors Qk belong 
to a discrete lattice. As the density increases the number 
of vectors Qk (and of domains around the Fermi sphere) 
may also increases[2] leading to a quasi-crystal which can¬ 
not be describe by a periodic model. 

Let us mention that in this paper we focus on the SDW 
states. These states are easier to compute leading to 
simpler formulas since the density of charge is constant. 
Equivalent results may be obtained[2l [3j for the charge 
density waves (CDW). 

In the following, we outline the main steps in the cal¬ 
culation of the SDW energies. First we introduce the de¬ 
formation of the Fermi surface generalizing Overhauser’s 
model and describe the SDW ansatz for the single parti¬ 
cle states. We then show how the optimal solution can 
be found by calculating the fixed point solution of a non¬ 
linear functional equation. The explicit results are then 
obtained by restricting to a one dimensional function and 
compared to the outcome of previous numerical simula¬ 
tions. 


Fermi gas energy of the truncated sphere. Let us 

call Afg the HF-energy of the Fermi gas where only plane 
wave states of wave vectors k inside the Fermi sphere 
of radius kp are occupied. Following Overhauser, in a 
first step the Fermi sphere is deformed into a volume 
J- as shown in Figjl] and its energy increase is denoted 
= EpQ — EpQ. Here, the subscript FG is used 
to point out that the many-body state is a Slater de¬ 
terminant of plane wave states inside the corresponding 
Fermi surface. Using kp as unit of wave vectors, the 
sphere in Figjl] has unit radius, and the deformation is 
characterised by a small parameter e approaching zero as 
Tg decreases. In order to keep the electron density con¬ 
stant, the deformed surface in the figure must be scaled 
byi? [U such that dk = R^ dk = 47r/3. The 

Fermi gas energy per particle in Hartree units (Ha = 


2 



FIG. 1: Shape of occupied zones for spiral SDW. J-^ (plain 
line) and (dashed line) are built from truncated sphere plus 
a small cylinder. A perturbed spin-up state is a superposition 
of a spin-up plane wave e'*“’ with a spin-down plane wave 
ghk+q)r Q^2{l-e + he). 



FIG. 2: AEpQ as a function of h for different values of e at 
Ts = 4. 


h^/{ma\)) is 


EpQ — 

Kfg = 
^FG = 


aKR^ 


Kfg- 


ayR^ 


Vfg 


dk 


I 


dkdk! 


1 


,1/3 


|k-k'|2 

1/3 ^ 


( 1 ) 

( 2 ) 

(3) 


with av = 3 ^ (x) ; aKiav = 271^ (^) « 37.9. 

The energy change, ATI^q, can be computed by direct 
integration (see Figj^ and gives at the leading order in 

em- 
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with a = 2(h — 4)^ -|- I and 7 = In ^ -|—Por small 

Ts, e is small and thus 7 is large. The minimum of AE^q 
with respect to /i is for /i = 4 — « 4 ^ and, for small 

Ts, AEpQ is four times smaller than in the Overhauser 
case {h = 0). 


The spin density waves. In a second step, the 
SDW are obtained by replacing a plane wave |k, f) by 
Ok |k,t) -f 6k |k -f Qk,i> (jokP + |6kP = 1 ), for k in JFf. 
Symmetrically, for k in J 7 , a plane wave |k,4,) is re¬ 
placed by the combination aj^ |k, 1 ) -f |k -|- Qk, t) with 
ajj. = a_k and = 6 _k. We choose Uk real and positive 
and in the following we assume that 6 k’s are also positive 
and that 6 k < av[T^. 

As in Fig[^ Qk is such that for k in J 7 , k + Qk does 
not belong to J 7 , and for k in J 7 , k-|-Qk does not belong 
to J 7 . The energy change is given by: 
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with kz = Q — kz (see Figj^. Using the linear symmetric 
operators T*: 


(r±/)(k) = ^dk'|^- 


k-k' 
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ik-k' 


/(k'), ( 8 ) 


Eq[5] rewrites: 

(2(«,62) + (T-a2,62)-(T+a6,a6)) (9) 

F s 


where {f,g) is the scalar product and 2 K{k) = 

2 av^ Q(Q — 2kz) > 0. In Eql^the difference between 1 
and R is negligible and in the lollowing we set i? = 1. 

Optimal solution. From the variations of Eq|^ with 
respect to 6 k, the optimal function 6 satisfies: 


26/c -I- 6 T {a? 


b^) = 



62 


r+a 6 


( 10 ) 


Setting ^ = ab and using 6 ^ < 4, — 6 ^ = 

Eq[^ rewrites ^ = J{0 where: 


AO 
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Thus the point is now to find the fixed points of the 
operator J. The Fermi gas (^ = 0) is a trivial fixed 
point. By definition 0 < ^ and from Eqllll we see 
that 0 < J(^) < 5 . We claim that starting witli^ = 1/2 
and iterating the process ^ —>■ J(^) leads to a non-trivial 
fixed point satisfying: 

akz=Q/2) = ^. 


( 12 ) 
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Indeed, by ([^ the kernels of are positive, thus 
are positivity preserving linear operators: if ^ then 
> r+^'; similarly < T-yJl/A - , 

and consequently: 

e > r ^ j{s,) > AO- (13) 

Thus starting with = 1 / 2 , we have J(Co) < Co and set¬ 
ting Cn = T(Cn-i), Cra is a decreasing sequence of positive 
functions and thus converges to a fixed point Coo- 

1-D approximation. Now we impose that (thus 
Ck) is non zero only in the cylinder C corresponding to 
the gray region of Fig{^ where it depends only on kz- 
C = {k : kl + ky < = I - {1 - e)^ 2e,0 < kz < Q}. 
As we shall see below, differs from zero only in the 
neighborhood of the top disk of (and its symmetric 
for J/). In any case, these restrictions always provides 
an upper bound for the energy of the SDW. 

First, for the second term of Eq[^ we have: 



FIG. 3: Renormalized energy AArg/e^ as a function of the 
density. The dashed-dotted line stands for the analytical solu¬ 
tion: AEVg/e^ = —0.115. A stand for the SDW simulations. 
Others symbols stand HF energies (see Fig.5 of [^). 


where the scalar product is now (/, g) = f{x)g{x), 

and become in terms of x: 


{A,T-b^) = {b^,T-A) = (&^^■l) - {b^,T-b^) (14) 


From Eq|^ T~1 = ujF(k) — ujF(k) where vjr is the po¬ 
tential induced by the truncated sphere. In the spherical 
case, the potential of the unit sphere is given by: 


?;(k) = 27r -|- TT 
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In this case, for k close to 1 (k and k are close and near 
the unit sphere and kz =2 — kz): t’(k) — u(k) « — 47 r(l — 
fc)log(l^) . This singular behavior is associated to the 
discontinuity of the density (in k-space). An analytic 
solution is provided for the truncated sphere [13) . This 
solution has the same behavior except that 1 — fc has to 
be replaced by the distance of k to the discontinuity of 
the density, i.e. the top disk of J/: 


/■r/r 
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Jo 

G{x) / dqdq' 

J n'^ ,q'^ <r'^ 


2<^2 r2x2-b (q - q')2 


( 20 ) 


=2 In 


1 + 1 


\x\u 
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In fact, for small r^, the term {T~b^,b^) may be ne¬ 
glected in Eqj^ [13]. In any case, since T~ is a positive 
operator, we get an upper bound for the energy and the 
variation of the resulting upper bound leads to ^ 
where J is now an operator on the positive functions on 
K+: 




1 T+C 

^ \J 7 r 2 {^x — 2x log(x))^ -b (T +^)2 


( 22 ) 


(T-l)(k)«-4^(Q/2-fc,)log(^^^/^2-^) (15) 


provided that |(5/2 — fc^j <C 1. For fc > 0, Eqj^is still 
valid [13] except in a small neighborhood of the edge of 
the top disk. In the following we neglect this effect and 
apply Eq[^also for fc > 0 . 

Using the scaled distance x = {Q/2 — kz)lr, 


2k+ T 1 = 27rr(7X — 2x log(x)), 
and integrating over q = (fc^,, ky), Eqj^ becomes: 






= 2^(7x- 2xlog(x),&2) - {T-b\b^) - {T+i,0 


(16) 
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As above, the fixed point of Eq[^can be easily found by 
iteration. 

Thereafter, for fixed the total energy variation 
A£’(e, fc) = AF;|Q(e, fc) -b is computed and op¬ 

timized with respect to e. For Cg = 3 about 20 iterations 
of the operator J are required and about 100 iterations 
for Tg = 0.01 (see A symbols in Figj^. In the next para¬ 
graph, we give a solution for ^ at small and deduce 
the scaling of AE from it. 

Analytic solution for small Cg. For small rg, 7 is 
large and Eqj^can be solved approximately m- 
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FIG. 4: Scaled parameter e/eo as a function of the density, eo 
is the value of Eq|28|for h = 1/2. The black circles stands for 
the present work, while other symbols stand for HF results[6]. 



FIG. 5: Influence of h on the energy gain verus rg: (a) energy 
and (•) e ratios. 


for 27r-\/2(7' — 7 ) = ^^(Tr^+d), leading to the asymptotic 
behavior of AE^- 




7exp(-^v^) (25) 


with C = 8e 

Now Eqj^and Eqj^provide the behavior of AE{e, h): 


AE — AE^^ + AEpQ 
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The minimum energy is at e = eo(l + 0(yT7)): 
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discretized Fermi sphere of 64^, 96^ and 128^ values of 
k. This corresponds to 32, 48 and 64 equally distributed 
values of k in the interval (0,1). For < 5 evidence for 
SDW ground states have been found [ 6 ]. In Figs and 
1 ^ we show the corresponding energy gain per number of 
SDW (2 for hexa up to 12 for bcc). The larger energy 
gain of the HF simulations for > 3 can be mostly at¬ 
tributed to a smoother and better optimized shape, com¬ 
pared to the simple cylinder used in the analytical SDW; 
other assumptions such as the 1 -dimensional approxima¬ 
tion decrease the energy further by a factor of 2 — 3. For 
Tg < 3, the discretization of the Fermi sphere becomes 
crucial, and even the simulations with 128^ k-points used 
in are insufficient to resolve the expected SDW am¬ 
plitudes leading to the standard Fermi gas ground state 
(e = 0 and = 0)- Direct numerical simulations of the 
SDW in this high density region will require a significant 
increase of k-points by several order of magnitudes. 


where 70 = uk/{ av’^T's). Eqj^ shows that at small Tg, 
AEr^/e^ goes to a constant. Figj^ shows the numerical 
results for the scaled energy at h = 1/2 (A symbols). 
This scaled energy is of order of -0.1 over a wide range 
of Cg. On the other hand, while Cq, Eqj^ varies over 
decades when rg decreases, Fig|^ shows that the ratio 
e/eo is a slowly varying function. The analytical result is 
supposed to be relevant for large 7 that is for ^ 1 . 
This can be verified on the figure: the next corrections 
in Eqs{^ and [M| behave as 

Influence of h. The dependency in h is through the 
parameter a, see Eq|^ Figj^ shows the effect of h on e 
and AE obtained numerically. For small rg, at /i = 1/2 
(thus a = 1 / 6 ), the energy is actually 16 times larger 
than in the Overhauser model (h = 0, a = 2/3). At 
larger Cg, this ratio can be significantly increased, e.g. 
it is about 200 for rg = 5; in this region we expect the 
energy gain to deeply rely on the precise shape of E. 

Comparison with HF simulations. In previous HF 
computations of the jelliumj^ |7j we have considered a 


Conclusion. Considering the ground state of the jel- 
lium in the Hartree-Fock approximation, we quantified 
the energy of the SDW suggested by Overhauser. Fur¬ 
thermore, we prove that a modification of the truncated 
Fermi sphere leads to an energy gain 16 to 200 larger 
than in the Overhauser model. 

Our results readily extends to a polarized model : in 
Eq 0 we have to take into account the direct potential 
which appears with a factor 1 /Q^ and thus is negligible 
at small rg (of order e^). 

In order to obtain the energy of jellium, the results 
of Figj^must be multiplied by the number of SDW. For 
simple periodic models considered in previous works, this 
factor varies from 2 (hexa) up to 12 (bcc). At very small 
Tg, the perturbation of the SDW is localized in tiny re¬ 
gions which do not interact, thus, one may suppose that 
we can have many of them distributed around the Fermi 
sphere giving rise to a quasi-periodic behavior of the den¬ 
sity. 
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